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ABSTRACT 

The fluctuation of the angular positions of reference extragalactic radio and optical 
sources under the influence of the irregular gravitational field of visible Galactic stars is 
considered. It is shown that these angular fluctuations range from a few up to hundreds 
of microarcseconds. This leads to a small rotation of the celestial reference frame. The 
nondiagonal coefficients of the rotation matrix are of the order of a microarcsecond. 
The temporal variation of these coefficients due to the proper motion of the foreground 
stars is of the order of one microsecond per 20 years. Therefore, the celestial reference 
frame can be considered inertial and homogeneous only to microarcsecond accuracy. 
Astrometric catalogues with microarcsecond accuracy will be unstable, and must be 
reestablished every 20 years. 
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INTRODUCTION 

Modern Very Long Baseline Interferometry (VLBI) obser- 
vations approach microarcsecond accuracy in the determi- 
nation of the positions of extragalactic radio sources. This 
accuracy will soon reach the fundamental limit of positional 
measurements associated with the curvature of space-time. 
Optical astrometric experiments currently under considera- 
tion (the GAIA project) will also be able to reach microarc- 
second accuracy. Here, we consider the effect of weak mi- 
crolensing on the rotation of the radio celestial reference 
frame based on radio quasars, though this effect is general, 
and will also be relevant, in principle, to the corresponding 
optical reference frame. 

Photons of a radio or optical source (a quasar, for in- 
stance) move along a path that is defined by the gravita- 
tional field of all the stars in our Galaxy. The stars in our 
Galaxy also move. As a result, their gravitational field is 
nonstationary. Thus, the paths of photons are not straight 
lines, but rather curves along nonstationary trajectories. 

Two photons that leave the quasar at different moments 
of time move along different paths. Therefore, the directions 
in space from which an observer detects these photons will 
be different. This causes fiuctuations in the position of the 
quasar. The aim of this paper is: 1) to calculate the value of 
variations of the angular positions of selected quasars due to 
this effect and 2) determine the influence of the gravitational 
fleld of visible Galactic stars on the stability of the celestial 
reference frame. 

The difference between the real and apparent angular 
positions is proportional to the mass of the body that bends 
the trajectory of the photon. The lower the mass of a body, 
the less the bending of the trajectory of a test photon. The 



main contribution to the effect, therefore, comes from the 
most massive and most dense populations of stars. 

Star-like objects can act on the propagation of light in 
two ways. The first is via microlensing, which forms two 
images (in place of the single true image) separated by 
an angular distance of the order of 10~^ arcsec. A close 
angular coincidence of the background quasar and deflect- 
ing body is necessary in this case. This microlensing effect 
can act as a space te lescope with extremely high resolution 
( Blair fc Sazhin 1993| ) . The microlensing effect was predicted 
by Paczinsky ( T986| ^ and later discovered by the MACHO, 
EROS, and OGLE groups (Alcock et al. 1993; Aubourg et 
al. 1993; Udalski et al. 1994). 

As an ordinary lens that has an infinite radius of grav- 
itational action, a foreground object can change the trajec- 
tory of photons from a background object even if the mi- 
crolensing effect is absent. Therefore, the second way star- 
like objects can act on the propagation of light is by produc- 
ing a small deviation of a photon trajectory from a straight 
line when the defiecting body is relatively far from the un- 
perturbed star-observer trajectory. A body with mass M 
and impact parameter p relative to the unperturbed trajec- 
tory bends a light beam by an angular amount 2M/p. The 
observer measures an angular deviation M/p from the un- 
perturbed position of the quasar. This idea was considered 
in Sazhin (1996) in connection with the limit of astromet- 
ric positional accuracy and more theoretically in Zhdanov 
(1995). Let us consider a population with a local density (in 
the vicinity of the Sun) of the order of n = 0.1 pc~'^ with 
masses for individual objects of the order of 0.1 M©. If an 
extragalactic object is observed, one can expect that at least 
one defiecting object will have an impact parameter of the 
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order of 



p = {-kRh ny 



Here, Rh is the radius of the Galactic halo. In this case, the 
observer expects to see an angular deviation of the appar- 
ent position of the order of 1 /iarcsec = /las from the true 
position of this extragalactic object. 

Only the second propagation effect noted above will be 
considered here. The first effect (microlensing) is a very rare 
event: it is necessary to observe one or two million stars in 
order to detect one microlensing effect. On the contrary, the 
second effect will act to some extent on every celestial source 
of light. One can expect a deviation of the order of 1 pas 
for every extragalactic source. Thus, the angular range of 
interest to us is from 1 pas to 1000 pas, which is well below 
the level of the microlensing effect, but is measurable using 



moderr VLBI techniques. The efi'ect of weak microlensing 



can already be probed at radio wavelengths due to the high 
accuracy of VLBI observations. 



WEAK MICROLENSING 

Let us consider a situation in which there are an extragalac- 
tic source (source of a test photon S), a massive body B 
with mass Mb (source of a gravitational field that bends 
the trajectory of the test photon), and an observer who de- 
tects the test photon (bent by the gravitational field of B). 
The direction of the incident photon does not coincide with 
the straight line connecting the extragalactic source and the 
observer. Therefore, the apparent celestial position of the 
extragalactic source does not coincide with its true position. 

The effect is defined by the Einstein cone, which has 
the value 

2 ^ 4GMb Lsb 

c2 Lob{Lsb + Lob) 

where Lqb and Lsb are the distances from light source S 
to the body B and from the observer O to the body B, 
respectively. In the microlensing effect, two images appear, 
which have angular distances from the body B: 



and 
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where ip is the angular distance from B to the true position 
of the quasar, and ip\ and ip2 are the angular distances of 
the two images relative to B. The brightness of the second 
image ip2 is inversely proportional to the fourth power of 
ifi, and in the case ip » ips, the observer cannot see the 
second image due to its weakness. In this case, the first im- 
age has approximately the same brightness as the intrinsic 
brightness, without the microlensing effect. The separation 
of the first image from the real position of the star is 
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The Galactic body B moves with some angular velocity 
Q. due to its peculiar motion in the Galaxy. In this case, one 
can write an equation for ip{t) as a function of time: 



where ipp is the angular impact parameter of the body rela- 
tive to the background quasar, and we take t = to be the 
moment of closest approach. It is easy to recalculate this 
equation from the modulus of the angular separation to the 
spherical coordinates a and S. 



VARIATIONS OF THE COORDINATES OF 
REFERENCE QUASARS 

The question of the variations of the coordinates of selected 
quasars is especially important because the 23rd General 
Assembly of the International Astronomical Union decided 
that, as of 1 January, 1998, the lAU celestial reference sys- 
tem is the International Celestial Reference System, in re- 



placement of the J2000 system realized by the FK5 (Fricke 
et al. 1988). The ICRS is realized by the International Ce- 



lestial Reference Frame (ICRF), defined by the equatorial 
coordinates of a set of selected extragalactic compact ra- 
dio sources determined using VLBI. The ICRF consists of 
a catalogue of eq uatorial coordinates of these radio sources 
( [McCarthy 1996^ . 

To test the effect of gravitational refraction (or mi- 
crolensing effect), we used the catalogue published in the 
1995 International Earth Rot ation Servic e (lERS) Annual 
Report (RSC(IERS)95 C02) ( [lERS 199£| ). This includes a 
total of 607 objects spread over the sky from declination 
—85° to -1-85°. The uncertainties in the coordinates are from 
50 -T- 2000 pas. There are 236 primary objects that are the 
most compact and best observed; 321 secondary compact 
sources that may have very precise coordinates in the future, 
when more observations are accumulated; and 50 sources 
that are complementary objects observed for optical frame 
ties or other objectives. 

The stability of the frame is based on the assumption 
that the sources have no proper motion. The hypothesis that 
the sources are fixed is used in performing the coordinate 
transformation between the celestial and terrestrial refer- 
ence systems. In reality, some of these sources have signifi- 
cant structure on mas scales. Changes in the source's bright- 
ness distribution can shift the effective brightness centroid 
of the source and thus its coordinates. A small number of ob- 



jccts from the reference frame exhibit such changes (Jacobs 
et al. 1993), and they can reach a few tenths of mas. 



Gravitational refraction leads to changes in the coordi- 
nates of all sources with time, and to a small rotation of the 
frame. Here, we do not calculate the effect of gravitational 
refraction from dark defiecting bodies and faint stars, but 
only from visible stars. The infiuence of dark bodies on the 
stochastic motion o f apparent ex tragalactic source positions 
was considered in (sazhin 1996), and is of the order of 1 
pas. The rotation of the celestial reference frame due to the 
influence of dark bodies and fainter stars will be considered 
in a later study. 

We used the Guide Star Catalogue (Lasker et al. 1990; 
Russel et al. 19 90; Jenkner et al. 1990 ) and the HIPPAR- 
COS catalogue (Ferryman et al. 1997) to find stars whose 



angular distances from the extragalactic reference sources 
do not exceed 2'. A total of 170 quasars have a total of 313 
nearby stars. 73 quasars have 86 neighboring stars in a circle 
with radius 1'. For most of these stars, we know neither their 
mass nor their distance. Two bright stars were found in the 
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HIPPARCOS catalogue. The proper motion, trigonometri- 
cal parallax, and spectral type of both stars are known, so 
that it is possible to determine their masses and distances. 

To estimate the rotation of the true ICRF relative to the 
observed ICRF, we used only the primary reference sources. 
27 sources have nearby stars, and 5 have two stars within 1'. 
The rotation vector = (Si, 02, S3) is defined by the matrix 




(1) 



where ri and ri' are the radius vector of the ith source be- 
fore and after the deflection of the light from this source. 
We estimated the di stance to th e star Lo b using the mass- 
luminosity relation (Allen 1973): 



where La , Lq are the luminosities and Ala , Mq the masses 
of the star and the Sun, respectively. The mass of each star 
was determined in the range from Mq to 30Mq using the 



Salpeter mass function dn/dM ~ 
can be found from the equation 



M- 



. The distance Lqb 



IgLoB + G.2AL0B = 0.2(m + 0.26) + 1.751g 



Ma 

Mq 



The apparent magnitudes m of the stars lie in the range 
from 11.3 — 15.8. The distances vary from 100 to 5500 pc. 
We assume here that the total absorption A in the vici nity 
of the Sun is ^ = 1.9"^ kpc'^ ( [Kaplan fc Pikelner 1979| ). 

We find that = (-0.08 ± 0.02, -0.47 ± 0.03, -0.36 ± 
0.03) pas for various realizations of the Salpeter mass dis- 
tribution (the variances of correspond to the ranges for 
various mass distributions). The change of with time was 
estimated using a uniform distribution of possible proper 
motions of the stars in the range — 100-^ 100 /^as/year. The 
value of the time derivative (or the angular rotation rate 
of the ICRF) is of the order of 1 ^as per 20 years, but can 
reach 1 pas per 1 year. The number of primary reference 
sources in the catalogue is 236. The resulting value of © 
is defined by the ~ 10% of the sources that have nearby 
stars. The stochastic process of ray deflection has an un- 
known distribution. If we increase the number of primary 
reference sources, the value of will be decreased, but the 
law by which diminishes is unknown. Since this law has a 
nongaussian distribution, the diminishing of will not be 
simply proportional to the inverse square root of the size of 
the catalogue. 

The largest effects of gravitational refraction were dis- 
played by several of the extragalactic sources (Table 1). The 
first row of Table 1 contains the name of the source, its 
right ascension and declination, and their uncertainties as 
reported by the lERS. The second row shows the magnitude 
and coordinates of the nearby star. In the third row, the an- 
gular distance (in arcsec) and corrections for the quasar's 
coordinates (in pas) are shown. These coordinate correc- 
tions were calculated assuming that the mass and distance 
of each star are Ma = SM© and Lqb = 500 pc. 

The exact effect of gravitational refraction can be cal- 
culated for the two stars from the HIPPARCOS catalogue 
(Table 2). The change of the stars' coordinates due to their 
proper motion leads to a change of the angular distance to 



n 
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Aa (mas) 

Figure 1. The variation of the coordinates of the quasar 0459-753 
due to the proper motion of the star 23106. 



the nearby quasar. As a result, a secular variation of the co- 
ordinates of the quasar will be observed. Figure 1 shows the 
change of the coordinates of the quasar 0459-753 over 100 
years, from 1952 to 2050. During the period from 1986-1996 
(when the observations were carried out) the right ascension 
and declination of the quasar changed by up to 7.5 mas and 
1.5 mas, respectively. The first star (with number 23106) has 
spectr al type Kl IIIp. It is a star with mass approximately 
4A/0 ( [Allen 1973| ). The parallax of this star is 3.43 ± 0.61 



mas and its proper motion is p^ = —4.00 ± 0.57 mas/year. 
Pi = —2.55 ± 0.78 mas/year. 

We suggest that the large uncertainties of the coordi- 
nates reported by lERS (Table 2) could be connected with 
the apparent motion of the quasar due to gravitational re- 
fraction. The calculated values of Aq and A5 are less than 
the observed values. The path of the quasar (Fig. 1) depends 
on the star's parameters, including its coordinates and mass. 
If the mass is underestimated, the true values of Aq and A(5 
should be larger. In the new version of the ICRF (avail- 
able at lERS under the label RSC(WGRF)95 R 01), the 
quasar 0459-753 is omitted (perhaps due to large coordinate 
uncertainties due to the unidentified effect of gravitational 
refraction) . 

Figure 2 shows the change of the coordinates of the 
quasar 1213-172 over 100 years. The second star (with num- 
ber 59803) has spectral type j? 8III. It is a normal giant 
with mass approximately 10 A/0 ( Allen 1973| ). The parallax 
of this star is 19.78 ± 0.81 mas and its proper motion is 
Pa = —159.58 ±0.66 mas/year, pa = 22.31 ±0.54 mas/year. 
The effect of refraction is smaller in this case, though the 
star is very close to the Sun (50 pc). 
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Table 1. Coordinates of quasars and nearby stars and the effect of gravitational refraction. Aa, A(5 are measured in /nas 



Source/Magnitude/ 
Distance (") 



Right ascension 



Declination 

1 55 



Uncertainty 



0007+106 
14.7 
0.358 



10 
10 



31.005871 
31.01 

Aa = 39.8 



10 
10 



58 
58 



29.50408 
29.9 

A(5 = 219.6 



0.000018 0.00042 



0111+021 
13.8 
0.509 



13 
13 



43.144949 
43.13 

Aa = -82.4 



22 
22 



17.31639 
17.8 

A5 = 133.7 



0.000014 0.00038 



0735+178 
4.3 

0.284 



38 
38 



7.393743 
7.37 

Aa = -293.6 



17 
17 



42 

42 



18.99868 
19.0 

A<5 = -30.5 



0.000003 0.00005 



0912+297 
15.3 
0.496 



15 52.401619 
15 52.40 

Aa = 18.1 



29 33 
29 33 



24.04293 
23.5 

A5 = -160.5 



0.000017 0.00034 



1101+384 
12.8 
0.054 



11 
11 



27.313911 
27.31 

Aa = 776.7 



38 
38 



12 
12 



31.79962 
31.8 

Ai5 = 1329.0 



0.000026 0.00038 



1302-102 
14.8 
0.812 



13 
13 



33.015008 
32.98 

Aa = -62.6 



-10 
-10 



33 
33 



19.42722 
20.1 

A5 = -77.0 



0.000018 0.00021 



1514^241 15 17 41.813132 -24 22 19.47552 

13.3 15 17 41.84 -24 22 19.8 

0.536 Aa = 120.8 A5 = -100.9 



0.000019 0.00031 



Table 2. The coordinates of the quasars and HIPPARCOS stars 



Source/Number 



Right ascension 
h m s 



Declination 



Uncertainty 

s " 



0459-753 
23106 



58 
58 



17.945614 
17.95 



-75 
-75 



16 
16 



37.95439 
38.0 



0.000923 
0.000037 



0.00315 
0.00060 



1213-172 
59803 



12 
12 



15 
15 



46.751743 
48.47 



-17 
-17 



31 
32 



45.40314 
31.1 



0.000043 
0.000042 



0.00041 
0.00049 



CONCLUSION 

There are two important conclusions from our investigation. 
First, we have shown that variations of the angular posi- 
tions of reference quasars range from a few microarcseconds 
up to hundreds of microarcseconds. The effect of weak mi- 
crolensing leads to a small rotation of the celestial reference 
frame. The nondiagonal coefficients of the rotation matrix 
are of the order of microarcsecond. The temporal variation 
of these coefficients due to stars' proper motions is of the or- 
der of one microarcsecond per 20 years. Because the proper 
motions of the majority of Galactic stars are unknown, the 
extragalactic reference frame can be considered inertial and 
homogeneous only to microarcsecond accuracy. The effect 
considered here should be investigated using current VLBI 
catalogues; however, this effect is general, and astromctric 
optical catalogues with microarcsecond accuracy will also be 
unstable, and must be reestablished about every 20 years. It 
is possible that one example of weak microlensing has been 



found: the source 0459-753, which has a nearby HIPPAR- 
COS star and shows large uncertainties in its coordinates. 
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Figure 2. The variation of the coordinates of the quasar 1213-172 
due to the proper motion of the star 59803. 
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